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Abstract. We present results for screening masses of mesons built from light and strange quarks in the 
temperature range of approximately between 140 MeV to 800 MeV. The lattice computations were per- 
formed with 2 + 1 dynamical light and strange flavors of improved (p4) staggered fermions along a line 
of constant physics defined by a pion mass of about 220 MeV and a kaon mass of 500 MeV. The lattices 
had temporal extents N T =4,6 and 8 and aspect ratios of N S /N T > 4. At least up to a temperature of 
140 MeV the pseudo-scalar screening mass remains almost equal to the corresponding zero temperature 
pseudo-scalar (pole) mass. At temperatures around 37^ {T c being the transition temperature) the contin- 
uum extrapolated pseudo-scalar screening mass approaches very close to the free continuum result of 2-kT 
from below. On the other hand, at high temperatures the vector screening mass turns out to be larger 
than the free continuum value of 2nT. The pseudo-scalar and the vector screening masses do not become 
degenerate even for a temperature as high as 4T C . Using these mesonic spatial correlation functions we have 
also investigated the restoration of chiral symmetry and the effective restoration of the axial symmetry. 
We have found that the vector and the axial-vector screening correlators become degenerate, indicating 
chiral symmetry restoration, at a temperature which is consistent with the QCD transition temperature 
obtained in previous studies. On the other hand, the pseudo-scalar and the scalar screening correlators 
become degenerate only at temperatures larger than 1.3T C , indicating that the effective restoration of the 
axial symmetry takes place at a temperature larger than the QCD transition temperature. 
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1 Introduction 

Hadron correlations at finite temperature have been advo- 
cated since long as a means to learn about non-perturbative 
properties of the quark gluon plasma from lattice simula- 
tions [1] . Studies of the in- medium properties of hadronic 
excitations provide information about some important length 
scales at high temperatures and give an idea about the 
relevant degrees of freedom in the plasma and their pos- 
sible physical effects. Furthermore, these studies also illu- 
minate aspects of the chiral and the effective anomalous 
J/a(1) symmetry restorations in QCD. 

Most lattice analyses of hadronic excitations have con- 
centrated on spatial correlation functions. Different from 



correlators in the temporal direction the physical extent 
of which is limited by the inverse temperature, a spatial 
correlation can be studied, in principle, at arbitrarily large 
distances facilitating the isolation of the lowest excitation. 
The exponential decay of these spatial correlators defines 
the so-called screening masses [1]. 

Physically, the inverse of a screening mass can immedi- 
ately be interpreted as the spatial distance beyond which 
the effects of putting a test hadron in the medium are 
effectively screened. The screening masses (M) need not 
be identical to the masses (m) defined via the exponen- 
tial decay of temporal correlators, but they are related to 
the same spectral function. Thus, amongst other purposes, 
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they provide a test bed for model building of the spectral 
function. 

Mesonic spatial correlation functions can be very use- 
ful tools to investigate the chiral and effective E/a(1) sym- 
metry restorations in QCD. Symmetry restorations de- 
mand that certain spatial (and also the temporal) correla- 
tion functions become identical at all distances. Of course, 
this also means that the screening masses of these mesons 
become degenerate as a symmetry gets restored. However, 
when one is working with a limited spatial extent, i.e. 
with relatively smaller volumes, extraction of screening 
masses from the large distance behavior of the correlators 
becomes difficult. Even in such cases one can obtain indi- 
cations of symmetry restorations by studying the degen- 
eracies between certain mesonic spatial correlation func- 
tions themselves at relatively shorter distances. Moreover, 
such symmetry restorations can be studied even without 
the use of the computationally demanding quark-line dis- 
connected diagrams. As for example, restoration of chiral 
symmetry can be studied through the degeneracy of the 
spatial correlation functions belonging to the vector (p) 
and the axial- vector (ai) channels. It is known that 
in the limit of mass-less light quarks the 2-point pseudo- 
scalar and scalar correlation functions of the 2 + 1 flavor 
theory are also sensitive to topologically nontrivial config- 
urations similar to the case of a 2 flavor theory. Hence the 
effective restoration of Ua(1) symmetry at high temper- 
ature can also be studied using the degeneracy between 
spatial correlation functions in the pseudo-scalar (7r) and 
the iso-triplet scalar (do) channels. 

Accurate non-perturbative determination of mesonic 
screening masses also plays an important role in testing 
the applicability of the dimensionally reduced perturba- 
tion theory [3]. Within the framework of dimensional re- 
duction at a particular temperature (T) one integrates 
out all the non-zero bosonic modes and all the fermionic 
modes of the 3+1 dimensional theory in order to obtain an 
effective 3 dimensional theory, for distance scales RT 3> 1 , 
consisting of static gauge fields coupled to adjoint scalars. 
Such a procedure is justified only when the spatial cor- 
relation lengths associated with all the fermionic modes 
in the original 3 + 1 dimensional theory are smaller than 
the largest spatial correlation length associated with the 
bosonic modes. Thus the applicability of dimensional re- 
duction, at a certain temperature and distance scale, cru- 
cially depends on the assumption that the screening masses 
of all the objects made out of fermionic modes of the 
original 3 + 1 theory are larger than the lowest screening 
mass associated with objects consisting only of the origi- 
nal gauge fields. In fact, this issue was addressed in Rcf. 
[1] in 4 flavor QCD using the standard staggered action on 
lattice with temporal extent N T = 4. This study suggests 
that the pseudo-scalar quark-anti-quark ('pion') correla- 
tor gives the largest spatial correlation length (i.e. small- 
est screening mass) for temperatures as high as T = 3T C , 
T c being the chiral transition temperature, indicating the 
breakdown of dimensional reduction below that temper- 
ature. On the other hand, another recent study [5] in 
2 + 1 flavor QCD using an improved staggered fcrmion ac- 



tion found that the non-perturbative values of the spatial 
string tension are in good agreement with the prediction of 
dimensionally reduced perturbation theory [5] right down 
to T = 1.5T C . Thus in order to clarify this issue related to 
the applicability of dimensional reduction it is necessary 
to have a detailed study of the meson screening masses 
in 2 + 1 flavor QCD using an improved staggered fcrmion 
action and small lattice spacings. 

In addition to the test of the applicability of the di- 
mensional reduction framework, improved determinations 
of mesonic screening masses are also important for check- 
ing the accuracy of the prediction of re-summed pertur- 
bation theories. In mass-less free field theory the meson 
screening masses acquire a value of 2irT [7] independent 
of the J PC structure of the interpolating meson operators. 
The leading order perturbative corrections of 0(g 2 T) to 
this free field value have been calculated using both the di- 
mensional reduction framework [8] and the hard thermal 
loop framework [5J. It was found that the leading order 
perturbative correction is positive and equal for all the 
meson channels. Thus for large temperatures the pertur- 
bative results approach the free field value from above in 
the same way independent of the meson channel. 

On the lattice, screening masses have been studied in 
the quenched approximation and with dynamical quarks, 
both in Wilson as well as staggered type fermion dis- 
cretizations (for a review and references see e.g. [lOj and 
references therein). Regarding the general pattern, little 
differences between quenched and full QCD studies have 
been observed. Below the transition region differences be- 
tween masses and screening masses were found to be small, 
independent of the lattice fermion action. Above the tran- 
sition mesonic screening masses were found to approach 
the zero quark mass continuum free field limit of 2nT from 
below. In the temperature range of 2 to 4 times the transi- 
tion temperature vector screening masses have been closer 
to this value while pseudo-scalar ones were found to be 
typically 10 % smaller than the vectors for Wilson quarks 
and somewhat lower values were reported for the standard 
staggered action. This is in contrast to the perturbative 
predictions as mentioned above. 

However, most of the studies were performed at fixed 
lattice spacing and fixed volume (V). More recently, the 
lattice spacing dependence of screening masses based on 
the standard staggered action was analyzed in [11]. The 
results show that discretization effects are significant for 
this choice of action and the pseudo-scalar screening mass 
at a fixed value of TV 1 / 3 approaches the infinite volume 
free lattice value from below. Preliminary results of a sys- 
tematic study of volume as well as lattice spacing depen- 
dence of the screening masses for quenched improved Wil- 
son fermions on large lattices have shown (T2"1[T5] that in 
this discretization scheme the volume effects increase the 
screening masses and discretization effects work in the op- 
posite direction. On the other hand, Ref. P3] concludes 
that finite volume as well as discretization effects are small 
for quenched overlap fermions, with the vector being close 
to 27rT and the scalar being 10% lower at T = 2T C . 
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Table 1. The list of meson operators studied in this work. 
States associated with the non-oscillating and the oscillating 
part of the screening correlators are designated by the identi- 
fiers NO and O respectively. Particle assignments of the corre- 
sponding states are given only for the u-d flavor combination. 



Here we present lattice results for meson screening 
masses in 2 + 1 flavor QCD using improved p4 [T5] stag- 
gered fermions with the zero temperature Goldstone pion 
mass of 220 MeV and the zero temperature kaon mass of 
500 MeV. The temperature ranges from about 140 MeV 
to maximally 820 MeV. Our study is performed at three 
different sets of the lattice spacings corresponding to N T — 
4, 6 and 8, which allows to address discretization effects. 

The paper is organized such that after presenting the 
numerical set-up in section [5] we briefly discuss discretiza- 
tion and finite volume effects as they arise, for our choice 
of action, in semi-analytical computations in the free case. 
This provides the background on which we present and 
discuss our results on screening masses in section 2] They 
have been determined from spatial correlators of operators 
built with degenerate up (u) and down (d) quarks as well 
as with strange (s) quarks with the same masses as in the 
generation of the configurations. Preliminary results have 
been presented earlier in [13lll6j . 



2 Numerical set-up 

In the staggered lattice formulation of QCD meson oper- 
ators can be defined [T7] as M = ip (r D ® r F ) tp, where 
the fermion field ip has four Dirac- and four flavor/taste- 
components and is defined on a coarse lattice with lat- 
tice spacing 2a. The matrices r D and r F are products of 
7-matrices and generate the spin- flavor/taste structure of 
the corresponding meson. Here we arc interested in local 
meson operators, for which r D = r F = r. In terms of 
staggered fermion fields {x) which are denned at the points 
x = (x, y, z, t) of a fine lattice with spacing a, these local 
meson operators can be written as M(x) = <p(x)\{x){x), 
where 4>{x) is a phase factor depending on the choice of 

r. 

In this work we study mesons built from ud, us and ss 
flavor combinations. For the flavor non-singlet combina- 
tions ud and us the correlation function consists solely of 
the quark-line connected part. However, for the flavor sin- 
glet ss flavor combination we only use the connected part 



of the correlation function and neglect the computation- 
ally demanding quark-line disconnected part. The con- 
nected part of the staggered meson screening correlator, 



projected to zero transverse momentum, p 
and to zero (boson) Matsubara frequency, uj, 
as 

c{z) = Y J <K*)({m£)*m£ 



= (Px,Py), 

is obtained 
(1) 



where Mq^ is the full staggered propagator from to 
x. Since a staggered fermion meson correlator, in general, 
contains two different mesons with opposite parity |18j the 
correlator is parametrized as 



C(z) 



Ano cosh 



M_ z 



Ns 
2 



(— l) z Ao cosh 







M + 


(-*)] 



(2) 



Different from e.g. [181119] we have chosen the phase fac- 
tor in Eq.© as <f>(x) = -{-l) x+y+T 4>{x). According to 
this convention M_ (M+) corresponds to the screening 
mass of the lightest negative (positive) parity stat^H and 
comes as the non-oscillating, NO (oscillating, O), part of 
the screening correlator. For instance, in this convention 
the Goldstone pion comes as the non-oscillating part of 
the screening correlator. Note that both amplitudes are 
positive PU, A N0 ,A > 0. 

For staggered fermions there are 8 possible local me- 
son operators [18] which are listed in Table [T] together 
with the corresponding phase factors 4>(x). We separate 
between the transverse (T) and longitudinal (£) polar- 
izations of vector and axial-vector correlators. At T = 
the rotational symmetry group of the lattice corresponds 
to the continuum 0(3). For spatial correlators, however, 
at non-zero temperature this symmetry breaks down to 
0(2) x Z(2) [311]. Thus, the longitudinal states are allowed 
to acquire screening masses different from the transverse 
ones. 

For the present study we have analyzed the gauge con- 
figurations generated by means of the RHMC algorithm 
[S] by the RBC-Bielefeld [MH and the HotQCD [53] 
collaboration using the p4 staggered action. These config- 
urations were generated for 2 degenerate dynamical light 
quarks and 1 heavier dynamical strange quark along the 
Line of Constant Physics (LCP). The LCP is obtained by 
tuning the bare quark masses such that at zero tempera- 
ture the (Goldstone) pion mass is approximately 220 
MeV and the kaon mass mx equals 500 MeV. As in [35] , 
relative scales have been obtained from the Sommer tq pa- 
rameter |24| . absolute ones from bottomonium splittings 
|25j resulting in ro = 0.469 fm in the continuum limit for 
physical quark masses. For more details about the simula- 
tions and scale settings confer Ref . [22] . The values of cou- 
pling constant and bare light quark mass in lattice units at 
which these configurations were generated are summarized 



1 For the screening states we keep, as usual, the J PC as- 
signments as they would emerge, for the same operators and 
correlation functions, in the zero temperature limit. 
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Table 2. Coupling constants /3 and light quark masses mi, 
temperatures T in MeV and the number of configurations of 
the given sizes on which screening masses were computed. The 
strange quarks have always been ten times heavier than the 
light ones, m s — 10m; . At the j3 values marked with a star we 
analyzed configurations which were generated in an early stage 
of the project at quark masses slightly off the LCP. 



in Table HJ The bare strange quark mass has always been 
10 times larger, m s = 10m;. The table also contains the 
temperature values and lists the number of configurations 
which were used for the present analysis. The configura- 
tions were separated by 10 trajectories. Autocorrelation 
times were estimated from the point-to-point spatial cor- 
relation functions at intermediate to large separations z 
from source and sink. The largest autocorrelations were 
observed in the lightest pseudo-scalar channel, M.2, and 
turned out to be of O(10) configurations, i.e. 0(100) tra- 
jectories, in the vicinity of the transition and dropping 
quickly away from it. 



3 Free case 

Lattice results are in general affected by discretization as 
well as by finite volume effects. At fixed temperature these 
are parametrized through the temporal extent of the lat- 
tice, aT = 1/N T , and the aspect ratio, TV 1/3 = N S /N T , 
respectively. Discretization and finite volume effects are 
readily computed semi-analytically in the non-interacting 
theory. While in this case both effects are expected to be 



larger than in the interacting theory [T!?1[T5] these compu- 
tations give at least hints about the sign and the order of 
magnitude of the lattice distortions. The free case is also 
the infinite temperature limit of the interacting theory. In 
this limit the spectral function underlying both, temporal 
as well as spatial correlation functions, solely consists of 
a cut in the complex frequency plane. The free case thus 
provides a test example for screening masses arising from 
a cut. 

In the free case at vanishing quark mass only states of 
one parity contribute to the correlation functions. In chan- 
nels M.2 and A45 to A47 these are states with negative 
parity while in the remaining channels the contributing 
states have positive parity. In the latter case an additional 
phase factor — (— l) z has been applied to C(z), see Eq.([2|) 
and Table [TJ such that for all channels the effective mass 
can be defined as 



aM eff (z) 



In 



C(z + 1) 
C(z) 



(3) 



With this sign convention the correlation functions of all 
channels are degenerate at zero quark mass for all values 
of N T and aspect ratios N S /N T > 4. At large z the spa- 
tial meson correlators start oscillating. These oscillations 
originate from a contribution ~ (— ) z ■ constant (in our 
phase convention) which vanishes in the thermodynamic 
limit N* —¥ oo. 



In Fig. 1(a) we show M e ff(z) as a function of the in- 
verse separation in units of temperature, l/(zT), for p4 
fcrmions at various values of the temporal lattice extent. 
As can be seen, for the improved p4 action there is practi- 
cally no N T dependence for zT > 1. On the other hand, for 
the unimproved, standard staggered action the discretiza- 
tion errors are significantly larger, as is shown in Fig. 1(b) 
for comparison. 

However, for both the p4 and the standard staggered 
action the effective mass from a point correlator is strongly 
dependent on z, and the value 2nT is obtained only in 
the limit of infinite z. This is in accord with the contin- 
uum, non-lattice calculation of the spatial pion correla- 
tion function [7] where this value is also only reached in 
the limit z — > oo of the continuum equivalent of Eq.Q. 
The z dependence of the effective screening mass is due 
to the cut in the spectral function. Because of contribu- 
tions from quarks at all momenta, p x ,Py and all Matsub- 
ara frequencies, uj n = nT(2n + 1) which die out only in 
the limit z — > oo the effective screening mass does not de- 
velop a plateau. This is strictly true for effective masses 
calculated from point-to-point correlation functions. Ex- 
tended sources, for instance wall sources, suppress contri- 
butions from higher quark momenta and thus may lead to 
a plateau like behavior in effective mass plots, at least in 
a limited z interval. Since infinite separations can only be 
achieved in the thermodynamic limit, Fig. [T] suggests that 
the screening masses determined from lattice simulations 
may have significant dependence on the maximum avail- 
able distance along which the correlation functions were 
measured. This can be considered as a finite volume ef- 
fect, although an indirect one as for aspect ratios larger 



5 




0.5 1 1/zT 1.5 2 0.5 1 1/zT 1.5 2 



(a) (b) 

Fig. 1. Dependence of the effective screening mass, Eq.©, on N T and on the separation zT , obtained from point correlators 
for a free theory with the p4 staggered action (a) and with standard staggered fermions (b). 



than 3 there is practically no volume dependence of the 
free results other than through the limited z range that 
is available. As a corollary, we note that the absence of a 
plateau in effective masses extracted from point correla- 
tors is a necessary condition to the verification of free field 
behavior in screening masses. 



4 Results 

The following results on screening masses of mesonic states 
built from light quarks have been obtained on dynamical 
Np = 2 + 1 gauge field configurations. The analysis was 
performed on lattices of size 16 3 x 4, 24 3 x 6 and 32 3 x 8, 
22 . 2]3] accompanied by some cross-checks on available high 
temperature lattices with larger aspect ratios at N T = 4 
and 6, (see Table [2]). On all those configurations we com- 
puted spatial correlations of meson operators (see TablcQ]) 
made out of degenerate u and d quarks. In addition we 
studied mesons with us and ss flavor content on the lat- 
tices with temporal extent N T = 4 and 6. In all cases the 
quark masses of the quark propagators have been chosen 
to be the same as in the quark determinant, i.e. in the 
generation of the configurations. On all lattices we com- 
puted point-to-point correlation functions. Subsequently, 
for most of the N T = 4 and 6 configurations, and some of 
the ones at N T = 8, (see the tables in the appendix) we 
also analyzed correlators from wall sources. 

In general, in a staggered fermion meson (screening) 
correlator a meson is accompanied by an opposite parity 
meson, see Eq. However, in the present study, in most 
channels we have seen indications for the presence of both 
of these mesons at most at the lowest temperatures. With 
increasing temperature, the contribution of one of the two 
mesons died out very fast. Hence, in most cases we found 
the signature of only one of the states. Specifically: 
- The M.1 channel is dominated by the oscillating positive 
parity scalar. A non-oscillating negative parity contribu- 
tion due to a pseudo-scalar state could not be isolated. 
We will denote the screening masses in this scalar channel 
as Msci keeping in mind that at least at zero tempera- 
ture the scalar channel contains two pseudo-scalar meson 



contributions which may become particularly important 
in the presence of taste violations [26] . 

- In the M2 channel we have not seen any oscillating 
contribution. This channel is dominated by the Goldstone 
pseudo-scalar at all temperatures. The screening masses 
in this channel will be denoted by Mpg. 

- The A43 and M.A. point correlators are noisy at low 
temperature but the signal-to-noise ratio improves greatly 
with rising T. Wall source correlators delivered suitable 
signals at all temperatures. In all fits we have found sta- 
ble signatures only of positive parity axial- vectors. The 
screening masses of the axial- vectors, denoted by Mav 
later, coming from these two channels are degenerate 

- In A46 and A47 we found only the negative parity vec- 
tor states. The signal-to- noise ratios were similar to the 
ones seen in A43 and M4. The screening masses of these 
states, My, coming from these two channels are found to 
be degenerate for all temperatures. 

- Channel A45 as well as A^8 were found to be very noisy 
in point correlators and reasonable signals could only be 
obtained from them at the highest temperatures. The wall 
correlators are at low temperature dominated by the axial- 
vector in Ai5 channel and the vector in Ai8 channel. At 
high temperature, independent of the source we found sig- 
natures of vector as well as axial-vector states in both 
channels. 

- These observations hold for all flavor combinations where, 
however, the signal-to-noise ratio is in general slightly bet- 
ter for the heavier combinations. 

- Contrary to the free case, in none of the correlation 
functions a contribution constant in z was observed. 

In all cases where acceptable signal-to-noise ratios were 
obtained within our statistics, the effective masses M e ff(z) 
from point-to-point correlators show reasonable plateaux 
in z within errors. This is particularly true at high temper- 
ature where we also analyzed lattices with aspect ratios 
larger than 4 to access larger zT values. There is a slight 
trend in the effective masses to decrease with increasing 
z, but in none of the channels the free field behavior as 



depicted in Fig. 1(a) could be verified. 

Subsequently we analyzed spatial meson correlators 
from wall sources built along the x, y, r directions. As 
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Fig. 2. (a) The pseudo-scalar screening masses in units of ro for iV T = 4, 6 and 8. The N T = 4 temperatures have been shifted 
by -8 MeV and the N T = 8 ones by +5 MeV, see text, (b) Flavor dependence of the pseudo-scalar screening masses at N T — 6. 
(c) Ratio of the pseudo-scalar screening masses at N T = 6 to the corresponding zero temperature masses determined at the 
same couplings along the line of constant physics with ~ 220 MeV and m,K — 500 MeV. 
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Fig. 3. (a) The scalar screening masses in units of ro for N T = 4, 6 and 8. Temperatures have been shifted as in Fig. 2(a) (b) 
Flavor dependence of the scalar screening masses at N T = 6. (c) Ratio of the scalar to the pseudo-scalar screening mass at high 
temperature. The data is from the N T — 6 and 8 lattices for the ud flavor combination and from N T — 6 for us and ss. 



pointed out in the previous section, effective masses from 
wall source correlators potentially pretend the presence of 
genuine bound states even in the free case because con- 
tributions from free quarks with higher momenta are sup- 
pressed by construction. However, the point-to-point cor- 
relation functions have already indicated the absence of 
a substantial contribution from a free quark cut at the 
temperatures investigated. 

The effective masses from the wall source correlators 
in general approach plateaux from below, as at zero tem- 
perature. For the scalar and the pseudo-scalar channels 
the error bars on the screening masses are of the same 
magnitude as for the point correlators. The wall source 
masses tend to be slightly lower than the point source 
ones but the differences are at the one a margin. Since 
the scalar and pseudo-scalar data set for point sources is 
most complete we will plot point source results in the fol- 
lowing, for wall source numbers we refer to the tables in 
the appendix. In the vector and axial- vector channels the 
signal-to-noise ratio is much better for the wall sources, in 
particular at low temperatures. At high temperatures the 
differences between wall and point source results are again 
small, with somewhat smaller error bars on the wall source 



results. Except for the N T = 8 lattices, we will therefore 
plot and discuss wall source results for those channels be- 
low. 

Independent of the source, the screening mass values 
which will be quoted in the following have been obtained 
from fits to the correlators with a minimum distance z m i„ 
from source and sink of about z m i n T ~ 1, except for 
the high temperature checks on lattices with aspect ratios 
larger than 4 where larger separations up to z m i n T ~ 2 
could be accessed. Errors were determined by jack-knifing 
the fits. 

For the purpose of showing our data as a function of 
a common temperature scale, throughout this paper the 
data for N T = 4 has been shifted in temperature by a value 
of AT = -8 MeV and the one for N T = 8 by AT = +5 
MeV to properly account for the N T dependence of the 
transition region [53] . Also throughout this paper we will 
use the chiral crossover temperature of T c ~ 196 MeV 
as our reference transition temperature. This corresponds 
to the chiral crossover temperature, determined from the 
peak of the disconnected chiral susceptibility, for 6 X 24 3 
lattices with ~ 220 MeV and m K ~ 500 MeV [23]. 
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In Fig. 2(a) we show our results for Mps, in the ud fla- 



vor combination, in units of ro for all N T values of 4, 6 and 
8 as a function of temperature. The pseudo-scalar screen- 
ing masses stay constant to a high degree of accuracy be- 
low the transition region and rise rapidly above. For in- 
stance, at a temperature of about 250 MeV the (lightest) 
pseudo-scalar mass has already reached a value of ~ 3.5T 
and at T ~ 400 MeV the value has increased to ~ 5.25T 
in temperature units. The rapid rise of Mps/T indicates 
that in this temperature range the scale for the pseudo- 
scalar screening mass is not set solely by the temperature. 



In Fig. 2(b) we compare the different flavor channels 
for N T = 6. Below the transition region, for all flavor 
combinations the pseudo-scalar screening masses stay ap- 
proximately constant in temperature and are markedly 
different from each other, in accord with the mass dif- 
ferences at zero temperature. Above the transition region 
the differences between the flavor channels become smaller 
with rising temperature and practically vanish at around 
300 MeV. 

The differences between zero temperature masses and 
finite temperature screening masses in the confined phase 
are inspected more closely in Fig. 2(c) As can be seen, 
the ratio of Mps to mps, the ordinary zero temperature 
Goldstone pion masses determined at the same couplings 
along the line of constant physics, starts differing from 
one for temperatures T > 170 MeV. This holds true for 
all the three flavor combinations ud, us and ss. However, 
in the range 170 MeV <T< 190 MeV, the differences are 
distinct but only of the order of 5 — 10% slightly depending 
on the quark mass. 

In Fig. 3(a) we show the temperature dependence of 
Msci normalized by ro, in the ud flavor combination for 
our lattices with temporal extents 4, 6 and 8. Note that 
Msc decreases with T at low temperature before it starts 
to ri se ra pidly at 200 MeV. This is seen more clearly in 
Fig. |3(b)| where we plot Msc from the N T = 6 lattices 
for all flavor combinations. Above 200 MeV the screening 
mass differences between the flavor channels persist up to 
about 325 MeV and vanish at higher temperature. 

In the presence of sufficiently strong taste violations 
and at small enough Goldstone boson masses, at least at 
zero temperature the scalar flavor non-singlct channel, do, 
may receive its dominant contribution from an unphysi- 
cal two-Goldstone boson threshold [55]. However, in our 
case of screening correlators we observe that at tempera- 
tures in the confined phase Msc as obtained from fits to 
the data is significantly larger than two times the Gold- 
stone boson screening mass. At the same time Msc is 
lower than twice the Kaon mass as well as the ir-q thresh- 
old where the rj mass was estimated with the Gell-Mann- 
Okubo formula rr£ = {&m 2 K — m^.)/3. This indicates that 
here indeed Msc denotes the screening mass of a genuine 
scalar. We then compare this scalar screening mass with 
the pseud o-scalar one at temperatures above 200 MeV in 
Fig." 



containing light quarks, i.e. the one with the light ud con- 
tent as well as the heavier us combination, agree with 
each other to a good degree of accuracy suggesting that 
those differences survive in the light quark chiral limit. 
The non-degeneracy of scalar and pseudo-scalar screen- 
ing mass still observed at temperatures up to at least 250 
MeV thus indicates that the anomalous E/a(1) symmetry 
is not effectively restored at this temperature. The non- 
zero mass splitting (Msc — Mps) for the ss sector for 
250 MeV < T < 325 MeV is most likely due to larger 
explicit Ua{X) breaking in presence of the larger strange 
quark masses. This is also supported by the fact that also 
in this sector this mass splitting vanishes for T > 325 
MeV which, as mentioned before, is the temperature range 
where the difference between the light and strange quark 
masses become unimportant. 

In Fig. 4(a) the results for My at different N T val- 



3(c) 



The figure highlights that the differences be- 
tween the two gradually disappear with rising tempera- 
ture, leading to degeneracy, but only above a temperature 
of 250 MeV. Note that the differences for the channels 



ues are shown. At temperatures below the transition the 
vector screening mass My in physical units is roughly con- 
stant in T , while for temperatures larger than about 250 
MeV its ratio with T, My/T, is approximately constant 
over the investigated temperature range. This indicates 
that different from the pseudo-scalar and scalar channels 
the temperature sets the scale for the vector screening 
masses already slightly above the transition. The ratio 
My/T is slightly larger than the continuum free result at 
infinite separation. However, a comparison with the free 
lattice result in the region around zT ~ 1 which typically 
is the lower boundary of the fit window shows that the 
data are about 20% lower than this number. Moreover, 
even at the highest temperatures available the effective 
masses from point correlators do not show the z depen- 
dence as in the free case. 

Fig. |4(b)| shows the flavor dependence of the vector 
screening masses as a function of temperature for N T = 6. 
At low temperatures the screening masses are ordered 
according to the strange quark content whereas already 
at a temperature of T ~ 225 MeV, the vector screen- 
ing masses are flavor independent, supporting that for the 
vector channel it is the temperature which sets the scale. 

We also inspect the degeneracy of Mps and My , which 
is expected at least in the limit of infinite temperature. In 
Fig. |4(c) we show the ratio Mps/My for the lightest fla- 
vor combination, noting that the ratio is independent of 
the quark flavor for temperatures larger than 300 MeV. 
Also for this ratio, the results from the N T = 4 lattices 
are somewhat lower than those from the other two lattice 
sets with N T = 6 and 8, which arc in good agreement with 
each other. Overall we find that Mps/My is statistically 
significantly lower than 1 even at temperatures T > 500 
MeV. Taking the N T = 8 result at T = 550 MeV at face 
value the difference from unity is about 9%. There is, how- 
ever, a slight trend that the ratio rises with temperature. 

Fig.[5]summarizes our results for the axial- vector screen- 
ing masses. Below the transition the screening masses are 
ordered according to their strangeness content but above 
the flavor dependence is weak, Fig. 5(b) Note that be- 
low 200 MeV the axial-vector screening mass is decreas- 
ing with rising temperature leading to a decrease in its 
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ratio to the vector screening mass, Fig. 5(c) This ratio 
becomes compatible with 1 around T m 200 MeV. Thus 
within errors axial-vector and vector screening masses be- 
come degenerate at this temperature and remain so at all 
larger temperatures. Although the quark masses are not 
at the chiral limit the restoration of chiral symmetry is 
thus also reflected in spectrum results. 

The longitudinal A45 and M8 channels in general have 
the largest error bars, also when wall sources are used. In 
the wall source correlators at temperatures below the tran- 
sition, M5 is dominated by an axial- vector while in M8 
only a vector could be isolated. Both states are, within er- 
rors, degenerate with the transverse ones. At high temper- 
ature, in wall as well as point correlators, where available, 
the vector contribution to A45 has the same screening 
mass as the axial- vector contribution to M8, M_(M5) = 
M+(At8). Both values are about 5% larger than the de- 
generate transverse screening masses, Mav or My. Simi- 
larly, the axial- vector contribution to A45 is equal to the 
vector contribution to M8, M+(M5) = M-(M8). These 
screening masses turn out to be roughly 20% larger than 
My- As already noted in Section [5] due to the breakdown 
of 0(3) rotational symmetry to 0(2) x Z(2) for spatial 
correlators at finite temperature the longitudinally 
polarized (axial-) vector screening states in A45 and A48 



do not need to be degenerate with the transverse polariza- 
tions of the V and AV channels. However, an explanation 
for the ordering observed here for the first time is not 
known to us. 



5 Summary and discussion 

Utilizing the gauge configurations generated by the RBC- 
Bielcfcld [51122] and the HotQCD [23] collaborations using 
the improved p4 fcrmion action along a line of constant 
physics, determined by ~ 220 MeV and raj- ~ 500 
MeV, we have investigated the screening masses of mesons 
in 2+ 1 flavor QCD from 8 different local meson operators 
listed in Table [1] In order to address discretization errors 
the analysis was carried out at three different values of 
the temporal lattice extent, N T = 4, 6 and 8. Although 
for most of our analyses we have used lattices with vol- 
umes TV 1 / 3 = 4, at high temperatures we have also ana- 
lyzed some lattices with larger volumes in order to access 
larger separations. Moreover, the quark mass dependence 
was studied by analyzing screening correlators for three 
different quark flavor combinations, viz. ud, us and ss. In 
order to gain some idea about the systematics we have 
computed correlation functions using both point and wall 
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Fig. 6. Effective screening masses in the pseudo-scalar chan- 
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respective screening masses obtained from cosh fits to the cor- 
relators. Except for the highest temperature, where the corre- 
lator was measured using wall source, all the correlators were 
obtained using point sources on N T — 6 lattices for the ud 
flavor combination. 



As discussed in Section [31 the presence of a cut in the 
spectral function will be manifested through a decreasing 
effective screening mass (coming from point source corre- 
lation function) as a function of increasing zT. Such a non- 
plateau-like behavior is a characteristic of the free theory. 
However, within our present level of accuracy and an avail- 
able range of 1 < zT < 2.5, wc could not identify a clear 
decrease of the effective screening masses with increasing 
zT and have rather observed plateau-like behaviors. As 
an illustration, in Fig. [S] wc show the effective screening 
masses in the pseudo-scalar channel as a function of zT 
for a large range of temperatures. As can be seen from 
this figure, we have not found any clear evidence of a non- 
platcau-likc behavior in the effective screening masses as 
seen in the free case (see Fig. [JJ). Moreover, the plateau- 
like structure in the effective screening masses are in good 
agreement with the values of the screening masses that 
have been obtained by fitting those correlators. 

In the low temperature region of 0.752^ < T < T c we 
have found that the PS screening mass remains almost 
constant (within < 10%) staying close to the correspond- 
ing zero temperature mass and then starts rising rapidly. 
In this low T regime also the V screening mass remains 
almost constant and then shows a rapid increase. On the 
other hand, the SC and the AV scalar screening masses 
show decreasing tendencies, attain minima around T c and 
then start rising again. 

In the intermediate temperature regime T c < T < 
1.5T C the spatial screening correlators show signals of chi- 
ral and effective J/a(1) symmetry restorations. As dis- 
cussed in Section [TJ chiral symmetry restoration is re- 
flected through the degeneracy of the V and AV correla- 
tion functions at all distances giving rise to degenerate V 
and AV screening masses. We found that, for N T = 6, My 
becomes approximately equal to May at a temperature 
T ~ 200 MeV. This temperature is approximately equal 



to the chiral transition temperature T c for these lattices 
identified via the peak of the disconnected chiral suscep- 
tibility. In order to gain some more insight into this de- 
generacy in Fig. 7(a) we show the correlation functions of 
the V and AV chlmTiels at T = 203 MeV, i.e. T ~ 1.1T C , 
for the N T = 6 lattice. Indeed we see that the signature of 
the chiral symmetry restoration is not only limited to the 
degeneracy of the screening masses, which are governed 
by the long distance behavior of the correlation functions, 
but also is reflected in the degeneracy between the cor- 
relation functions at very short distance scales zT < 1. 
Similarly, as also mentioned in Section [TJ the spatial cor- 
relation functions of the PS and the (iso-triplet) SC show 
effective restoration of J7a(1) for the 2 + 1 flavor theory. In 
the present work we have found that the PS and SC cor- 
relation functions become degenerate, and hence also the 
respective screening masses (see Fig. 3(c) I, for T > 1.3T C 
only. To illustrate this fact in detail we show comparisons 
of th e PS and SC screening corr elator s at T ~ 1.2T C in Fig. 
[7(b)] and at T ~ 1.3T C in Fig. ff(cj\ While at T ~ 1.2T C 
the PS and SC screening correlators start differing for dis- 
tances zT > 0.5 the correlators become degenerate at all 
distance scales, even for zT < 1, at T ~ 1.3T C . Moreover, 



as shown in Fig. 3(c) the magnitude of the mass splitting 
(M sc ~ M PS ) intEe temperature range T c < T < 1.3T C 
remains almost unchanged in going from ud flavor con- 
tent to us. This indicates that in this temperature regime 
the mass splitting (Msc — Mps), which is sensitive to the 
presence of topologically non-trivial configurations [5] , re- 
mains non-vanishing in the limit of mass-less light quarks. 
This conveys the point that the effective J7a(1) restoration 
in QCD docs not happen at the chiral transition tempera- 
ture T c , i.e. the temperature where the V and AV correla- 
tion functions become degenerate. This observation may 
have implications in relation to the possible existence of 
a first order chiral transition in the limit of two mass-less 
flavors [27]. 

As has been pointed out previously in Section [TJ at 
higher temperatures T > 1.5T C the mesonic screening 
masses play an important role in testing the applicabil- 
ity of the dimcnsionally reduced perturbation theory. At 
these temperatures we have observed that the screening 
masses in the PS (SC) channel are always lower than the 
V (AV) screening masses. In other words, the largest spa- 
tial correlation length belongs to the PS (SC) channel. In 
order for the dimensional reduction to work this PS spa- 
tial correlation length must be smaller than the smallest 
spatial correlation length associat ed wi th a glueball like 
object. As can be seen from Fig. |8(a) Mps/T reaches 
values greater than 4 at T ~ 1.5T C and greater than 5 
at T ~ 2T C for all three lattice spacings. Thus, already 
for T ~ 1.5T C the value of the 'pion' screening mass is 
larger than the smallest glueball screening mass estimated 
in Ref. [T| and larger than the Debye mass estimated in 
Rcf. [55]. Hence, the lightest fermionic mode is compar- 
atively heavy already at rather small temperatures such 
that dimensionally reduced QCD may work down to tem- 
peratures as low as about 1.5T C , as was demonstrated for 
a quantity like the spatial string tension [5,6 . This is in 
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contrast to the previous findings of Ref. [3] where it was 
concluded that the dimensional reduction framework may 
fail for T < 3T C . 



in Fig. 8(c) show extrapolations when all the three lattice 



In Figs. 8(a) and |8(b)] we show close-up views of how 



the PS and ffieV screening masses deviate from their con 
tinuum free value of 2itT in the high temperature region 
T > 1.5T C . In this temperature regime for the pseudo- 
scalar (and scalar) channel the ratio Mps/2irT not only 
rises rapidly towards 1 with increasing temperature, but 
also comes closer towards 1 with decreasing lattice spacing 
at a fixed temperature. For example, while on our coarsest 
lattice Mps deviates from 2nT by about 10% at a temper- 
ature as high as 4T C this deviation becomes ~ 5% for our 
finest lattice already at T ~ 3T C . On the other hand, the 
ratio Mv/^ttT shows a rather slow increase with increas- 
ing temperature. While this ratio remains ~ 2% below 1 at 
T ~ 4T C for our coarsest lattice spacing it becomes ~ 4% 
larger than 1 at T ~ 3T C for our finest lattice spacing. In 
view of this, in Fig. 8(c) we attempt linear in a 2 ~ 1/N% 
extrapolations of the PS and V screening masses at a tem- 
perature T ~ 3T C , i.e. the maximum temperature where 
we have data at all three lattice spacings. The dashed lines 



spacings are used during the extrapolations. Such extrap- 
olations give that the PS and the V screening masses, 
extrapolated to zero lattice spacing, are respectively 2-3% 
below and 5-6% above the free theory value of 2ttT. As 
even at this high temperature the lattice spacing corre- 
sponding to N T = 4 is still quite large (about 0.09 fm) it 
is not clear whether including these data in our zero lattice 
spacing extrapolations yields reliable results. In order to 
have some quantitative idea on the systematics we have 
also performed extrapolations linear in 1/N 2 using only 
N T = 6 and 8 data. These are indicated by the solid lines 
in Fig. 8(c) Such extrapolations give that the continuum 



PS screening mass at this temperature may even be closer 
to 2itT, smaller only by 1%, and the continuum V screen- 
ing mass is about 7% larger than 27tT. 

Perturbation theory predicts (SJEl that at high tem- 
peratures the screening masses of all the mesons are de- 
generate and larger than the free theory value of 2wT. 
Although we have found that the continuum extrapolated 
PS screening mass approaches the continuum free theory 
value of 2-kT very closely around T ~ 3T C and also, for 
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and the BlueGcne/L at the New York Center for Compu- 
tational Sciences (NYCCS). 

Appendix 

Tables [3] [4] and [5] summarize the meson screening masses 
for the scalar, pseudo-scalar, averaged transverse axial- 
vector and averaged transverse vector channels for the 
lightest ud flavor combination and for all the lattices that 
have been analyzed in this paper. The screening masses 
are given in lattice units. The data can be easily converted 
to either temperature or vacuum (tq) units by means of 
the included zero temperature results for ro/a from [55] . 

References 



the first time, found evidence that the V screening masses 
probably become larger than 2irT at high temperature, we 
have not found any evidence that Mps and My become 
degenerate up to temperatures at least as high as T ~ 4T C . 
In order to have a closer look at this non-degeneracy, in 
Fig. [5] we compare directly the PS and V correlation func- 
tions as function of the distance zT around T ~ AT C . It 
is clear from this comparison that the spatial correlators 
in the PS and V channels differ from each other at all 
distances including short ones of zT < 1. In contrast, in 
the (mass-less) free theory the PS and V correlation func- 
tions themselves are degenerate at all distances even on a 
finite lattice. Although this observation is subject to cut- 
off and finite volume effects which could, in principle, be 
very different between the free and the interacting theory, 
we take it as a strong indication that the non-degeneracy 
of M ps and My is probably not a lattice artifact and may 
arise due to the presence of spin dependent interactions, 
as suggested e.g. in [25], at the investigated temperatures. 

As mentioned earlier, one of the sources of systematic 
errors in our study is the limited value of the maximum ac- 
cessible distance zT. We have also found moderate lattice 
spacing dependence of screening masses at high tempera- 
ture. Hence, it is important to corroborate our findings by 
performing analyses on lattices with larger aspect ratios 
and finer lattice spacings. We hope to address this issue 
in the future. 
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ro/a 


aMsc 


aMps 


ciMav 


aM v 


16 a x 4 


point sources 


3.290 


1.823 


0.860(80) 


0.293(08) 






3.320 


1.908 


0.702(20) 


0.349(08) 






3.351 


2.069 


0.695(18) 


0.517(08) 


1.33(12) 


1.37(21) 


3.382 


2.230 


0.802(11) 


0.778(15) 


1.28(10) 


1.33(12) 


3.410 


2.503 


0.920(10) 


0.902(18) 


1.46(06) 


1.46(06) 


3.460 


2.890 


1.046(48) 


1.043(46) 


1.54(09) 


1.52(08) 


3.540 


3.687 


1.162(36) 


1.162(47) 


1.48(05) 


1.49(06) 


3.570 


4.009 


1.206(10) 


1.206(10) 


1.53(01) 


1.53(01) 


3.630 


4.651 


1.259(09) 


1.259(09) 


1.54(01) 


1.54(01) 


3.690 


5.201 


1.306(08) 


1.302(04) 


1.55(01) 


1.55(01) 


3.760 


6.050 


1.341(07) 


1.341(07) 


1.56(02) 


1.57(02) 


3.820 


6.752 


1.369(10) 


1.369(10) 


1.57(02) 


1.57(02) 


3.920 


7.590 


1.400(05) 


1.400(05) 


1.56(01) 


1.56(01) 


wall sources 


3.382 


2.230 


0.801(13) 


0.762(19) 


1.380(60) 


1.380(50) 


3.410 


2.503 


0.905(11) 


0.890(13) 


1.460(60) 


1.450(70) 


3.460 


2.890 


1.054(22) 


1.051(22) 


1.515(35) 


1.520(38) 


3.540 


3.687 


1.166(10) 


1.166(10) 


1.506(09) 


1.506(09) 


3.570 


4.009 


1.204(06) 


1.204(06) 


1.516(05) 


1.516(05) 


3.630 


4.651 


1.255(05) 


1.255(05) 


1.517(04) 


1.517(04) 


3.690 


5.201 


1.291(05) 


1.291(05) 


1.528(06) 


1.528(06) 


3.760 


6.050 


1.330(04) 


1.330(04) 


1.540(10) 


1.540(10) 


3.820 


6.752 


1.355(04) 


1.355(04) 


1.541(05) 


1.541(05) 


3.920 


7.590 


1.386(02) 


1.386(02) 


1.545(02) 


1.545(02) 


32 a x 4 


point sources 


3.920 


7.814 




1.395(2) 




1.560(1) 


wall sources 


3.920 


7.814 




1.392(1) 




1.554(2) 



Table 3. Screening masses from N T = 4 lattices. 
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p 


ro/a 


aMsc 


aMps 


clMav 


aM v 


24 a x 6 


point sources 


3.351 


2.069 


0.900(30) 


0.257(01) 






3.410 


2.503 


0.700(20) 


0.226(02) 






3.430 


2.647 


0.615(09) 


0.216(03) 






3.445 


2.813 


0.510(20) 


0.253(05) 






3.455 


2.856 


0.420(10) 


0.292(07) 






3.460 


2.890 


0.435(08) 


0.334(08) 






3.490 


3.223 


0.490(10) 


0.440(10) 






3.510 


3.426 


0.580(10) 


0.560(10) 






3.540 


3.687 


0.648(09) 


0.640(08) 


1.11(6) 


1.08(6) 


3.570 


4.009 


0.728(08) 


0.726(08) 


1.04(3) 


1.04(3) 


3.630 


4.651 


0.790(10) 


0.790(10) 


1.04(1) 


1.04(2) 


3.690 


5.201 


0.851(08) 


0.851(08) 


1.04(1) 


1.04(1) 


3.760 


6.050 


0.900(04) 


0.900(04) 


1.06(1) 


1.06(1) 


3.920 


7.590 


0.957(02) 


0.957(02) 


1.07(1) 


1.07(1) 


wall sources 


3.351 


2.069 


1.070(50) 


0.270(3) 


1.650(20) 


1.270(30) 


3.410 


2.503 


0.780(20) 


0.240(2) 


1.408(09) 


1.125(05) 


3.430 


2.647 


0.685(02) 


0.227(1) 


1.140(50) 


0.985(06) 


3.445 


2.813 


0.512(02) 


0.246(1) 


1.090(04) 


0.994(04) 


3.455 


2.856 


0.456(02) 


0.262(2) 


1.058(05) 


0.990(10) 


3.460 


2.890 


0.436(01) 


0.316(1) 


1.010(08) 


1.030(10) 


3.490 


3.223 


0.487(01) 


0.438(1) 


0.990(06) 


0.990(07) 


3.510 


3.426 


0.562(01) 


0.538(2) 


1.009(04) 


1.002(04) 


3.570 


4.009 


0.705(01) 


0.702(1) 


1.016(01) 


1.011(01) 


3.630 


4.651 


0.788(01) 


0.788(1) 


1.031(01) 


1.031(01) 


3.690 


5.201 


0.839(01) 


0.839(1) 


1.043(02) 


1.043(02) 


3.760 


6.050 






1.049(05) 


1.050(04) 


3.920 


7.590 


0.938(01) 


0.941(1) 


1.059(01) 


1.059(01) 


32 a x 6 


point sources 


3.820 


6.864 




0.916(3) 




1.059(4) 


3.920 


7.814 




0.948(2) 




1.061(4) 


4.000 


9.048 




0.965(3) 




1.063(4) 


4.080 


10.390 




0.980(3) 




1.064(2) 


wall sources 


3.820 


6.864 




0.911(2) 




1.047(2) 


3.920 


7.814 




0.946(1) 




1.059(1) 


4.000 


9.048 




0.961(1) 




1.062(2) 


4.080 


10.390 




0.978(2) 




1.066(2) 



Table 4. Screening masses from N T = 6 lattices. 
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ro /a 


aMsc 


aMps 


qMav 


aM v 


32 a x 8 


point sources 


3.430 


2.647 




0.201(01) 






3.500 


3.328 


0.463(174) 


0.175(01) 






3.570 


4.009 


0.366(010) 


0.311(14) 






3.585 


4.160 


0.382(010) 


0.351(10) 






3.760 


6.050 


0.615(008) 


0.614(08) 


0.777(12) 


0.776(12) 


3.820 


6.864 


0.666(006) 


0.666(06) 


0.807(07) 


0.807(07) 


3.920 


7.814 


0.702(003) 


0.702(03) 


0.802(04) 


0.802(04) 


4.000 


9.048 


0.721(001) 


0.722(02) 


0.803(02) 


0.805(03) 


4.080 


10.390 


0.740(003) 


0.740(03) 


0.813(03) 


0.813(03) 


4.200 


12.420 


0.753(006) 


0.753(05) 


0.831(12) 


0.816(09) 


wall sources 


4.000 




0.716(001) 


0.716(03) 


0.799(04) 


0.800(03) 


4.200 




0.755(002) 


0.755(01) 


0.832(12) 


0.818(02) 



Table 5. Screening masses from iV T = 8 lattices. 



